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Abstract. The past, present and future of cosmic microwave background (CMB) 
anisotropy research is discussed, with emphasis on the Boomerang and Maxima bal- 
' loon experiments. These data are combined with large scale structure (LSS) infor- 

t — , mation derived from local cluster abundances and galaxy clustering and high red- 

shift supernova (SN1) observations to explore the inflation-based cosmic structure for- 
mation paradigm. Here we primarily focus on a simplified inflation parameter set, 
{u>b, u>cdm, &tot, ^Qj wq, n s ,Tc,crg}. After marginalizing over the other cosmic and ex- 
, perimental variables, we find the current CMB+LSS+SN1 data gives Vt to t = 1.04±0.05, 

consistent with (non-baroque) inflation theory. Restricting to Cltot = 1, we find a nearly 
scale invariant spectrum, n s = 1.03 ± 0.07. The CDM density, uj c dm = 0.17 ± 0.02, is 
in the expected range, but the baryon density, uib = fibh 2 = 0.030 ± 0.004, is slightly 
larger than the current 0.019 ± 0.002 Big Bang Nucleosynthesis estimate. Substantial 
dark (unclustered) energy is inferred, £Iq « 0.68 ± 0.05, and CMB+LSS £1q values 
are compatible with the independent SN1 estimates. The dark energy equation of 
state, parameterized by a quintessence- field pressure-to-density ratio wq, is not well 
determined by CMB+LSS (w Q < -0.3 at 95% CL), but when combined with SN1 
' the resulting wq < —0.7 limit is quite consistent with the wq—— 1 cosmological con- 

stant case. Though forecasts of statistical errors on parameters for current and future 
experiments are rosy, rooting out systematic errors will define the true progress. 
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CMB ANALYSIS: PAST, PRESENT AND FUTURE 

The CMB is a nearly perfect blackbody of 2.725 ± 0.002 K [1], with a 3.372 ± 
0.007 mK dipole associated with the 300 km s -1 flow of the earth in the CMB, and 
a rich pattern of higher multipole anisotropies at tens of /iK arising from fluctua- 
tions at photon decoupling and later. Spectral distortions from the blackbody as- 
sociated with starbursting galaxies detected in the COBE FIRAS and DIRBE data 
are due to stellar and accretion disk radiation being downshifted into the infrared by 



dust then redshifted into the sub millimetre; they have energy about twice all that 
in optical light, about a tenth of a percent of that in the CMB. The spectrally well- 
defined Sunyaev-Zeldovich (SZ) distortion associated with Compton-upscattering 
of CMB photons from hot gas has not been observed with FIRAS, but only at high 
resolution along lines-of-sight through dozens of clusters — with very high signal- 
to-noise though. The FIRAS 95% CL upper limit of 6.0 x 10~ 5 of the energy in the 
CMB is compatible with the < 10~ 5 expected from clusters, groups and filaments 
in structure formation models, and places strong constraints on the allowed amount 
of earlier energy injection, e.g., ruling out mostly hydrodynamic models of LSS. 

Upper Limit Experiments from the 70s & 80s: The story of the experi- 
mental quest for anisotropics is a heroic one. 1 The original 1965 Penzias and Wilson 
discovery paper quoted angular anisotropies below 10%, but by the late sixties 10~ 3 
limits were reached, by Partridge and Wilkinson and by Conklin and Bracewell. 
As calculations of baryon-dominated adiabatic and isocurvature models improved 
in the 70s and early 80s, the theoretical expectation was that the experimental- 
ists just had to get to 10 -4 , as they did, e.g., Boynton and Partridge in 73. The 
only signal found was the dipole, hinted at by Conklin and Bracewell in 73, but 
found definitively in Berkeley and Princeton balloon experiments in the late 70s, 
along with upper limits on the quadrupole. Throughout the 1980s, the upper limits 
kept coming down, punctuated by a few experiments widely used by theorists to 
constrain models: the small angle 84 Uson and Wilkinson and 87 OVRO limits, 
the large angle 81 Melchiorri limit, early (87) limits from the large angle Tenerife 
experiment, the small angle RATAN-600 limits, the 7°-beam Relict-1 satellite limit 
of 87, and Lubin and Meinhold's 89 half-degree South Pole limit, marking a first 
assault on the peak. 

These upper limit experiments were highly useful, in particular to rule out adi- 
abatic baryon-dominated models. In the early 80s, dark matter dominated uni- 
verses lowered theoretical predictions by about an order of magnitude. In the 84 to 
mid-90s period, many groups developed codes to solve the perturbed Boltzmann- 
Einstein equations when dark matter was present. Armed with these pre-COBE 
computations, plus the LSS information of the time, a number of very interesting 
models fell victim to the data: scale invariant isocurvature cold dark matter models 
in 86, large regions of parameter space for isocurvature baryon models in 87, infla- 
tion models with radically broken scale invariance leading to enhanced power on 
large scales in 87-89, CDM models with a decaying (~ keV) neutrino if its lifetime 
was too long (> lOyr) in 87 and 91. Also in this period there were some limited 
constraints on "standard" CDM models, restricting Vl tot) Q B , and the amplitude 
parameter as- (cr| is a bandpower for density fluctuations on a scale associated 
with rare clusters of galaxies, 8h _1 Mpc, where h = Hq/(100 km s -1 Mpc -1 ).) 

Post-DMR Experiments: The now familiar motley pattern of anisotropies 
associated with 2 < £ < 20 multipoles at the 30fiK level revealed by COBE at 
7° resolution was shortly followed by detections, and a few upper limits (UL), at 
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higher t in 19 other ground-based (gb) or balloon-borne (bb) experiments — most 
with many fewer resolution elements than the 600 or so for COBE. Some predated 
in design and even data delivery the 1992 COBE announcement. Proceeding from 
the period we began analyzing them, we have the intermediate angle SP91 (gb), 
the large angle FIRS (bb), both with strong hints of detection before COBE, then, 
post-COBE, more Tenerife (gb), MAX (bb), MSAM (bb), white-dish (gb, UL), 
argo (bb), SP94 (gb), SK93-95 (gb), Python (gb), BAM (bb), CAT (gb), OVRO- 
22 (gb), SuZIE (gb, UL), QMAP (bb), VIPER (gb) and Python V (gb). A list 
valid to April 1999 with associated bandpowers is given in [10], and are referred 
here as 4.99 data. They showed evidence for a first peak [10], although it was not 
well localized. Within limited parameter sets, good constraints on n s , some on Q tot 
and Q\ could be given, when LSS was added. 

The Present, TOCO, BOOMERANG & MAXIMA: The picture dramat- 
ically improved this year, as results were announced first in summer 99 from the 
ground-based TOCO experiment in Chile [2], then in November 99 from the North 
American balloon test flight of Boomerang [3]. These two additions improved peak 
localization and gave evidence for Q tot ~ 1. Then in April 2000 results from the 
first CMB long duration balloon (LDB) flight, Boomerang [4], were announced, 
followed in May 2000 by results from the night flight of Maxima [5]. Boomerang's 
best resolution was 10', about 40 times better than that of COBE, with tens of 
thousands of resolution elements. Maxima had a similar resolution but covered an 
order of magnitude less sky. 

Boomerang carried a 1.2m telescope with 16 bolometers cooled to 300 mK in the 
focal plane aloft from McMurdo Bay in Antarctica in late December 1998, circled 
the Pole for 10.6 days and landed just 50 km from the launch site, only slightly 
damaged. In [4], maps at 90, 150 and 220 GHz showed the same spatial features 
and the intensities were shown to fall precisely on the CMB blackbody curve. The 
fourth frequency channel at 400 GHz is dust-dominated. Fig. 1 shows the 150 GHz 
map derived using only one of the 16 bolometers. Although Boomerang altogether 
probed 1800 square degrees, only the region in the rectangle covering 440 square 
degrees was used in the analysis described in [7,8] and this paper. Fig. 1 also shows 
the 124 square degree region of the sky (in the Northern Hemisphere) that Maxima- 
1 probed. Though Maxima was not an LDB, it did so well because its bolometers 
were cooled even more than Boomerang's, to 100 mK, leading to higher sensitivity 
per unit observing time, it had a star camera so the pointing was well determined, 
and, further, all frequency channels were used in creating its map. 

Primary CMB Processes and Soundwave Maps at Decoupling: Both 
Boomerang and Maxima were designed to measure the primary anisotropies of 
the CMB, those which can be calculated using linear perturbation theory. What 
we see in Fig. 1 are, basically, two images of soundwave patterns that existed 
about 300,000 years after the Big Bang, when the photons were freed from the 
plasma. The visually evident structure on degree scales is even more apparent 
in the power spectra of the Fourier transform of the maps, which show a domi- 
nant (first acoustic) peak, a less prominent (or non-existent) second one, and the 



possible hint of a third one from Maxima. Fig. 1 also shows that the quite het- 
erogeneous 4.99+TOCO+Boomerang-NA mix of CMB data is very consistent with 
what Boomerang-LDB and Maxima show. 

The images are actually a projected mixture of dominant and subdominant phys- 
ical processes through the photon decoupling "surface", a fuzzy wall at redshift 
z r ~ 1100, when the Universe passed from optically thick to thin to Thomson 
scattering over a comoving distance ~ 10 h -1 Mpc. Prior to this, acoustic wave 
patterns in the tightly-coupled photon-baryon fluid on scales below the comoving 
"sound crossing distance" at decoupling, < 100 h -1 Mpc (i.e., < 100 kpc physical), 
were viscously damped, strongly so on scales below the ~ 10 h -1 Mpc thickness 
over which decoupling occurred. After, photons freely-streamed along geodesies to 
us, mapping (through the angular diameter distance relation) the post-decoupling 
spatial structures in the temperature to the angular patterns we observe now as the 
•primary CMB anisotropies. The maps are images projected through the fuzzy de- 
coupling surface of the acoustic waves (photon bunching), the electron flow (Doppler 
effect) and the gravitational potential peaks and troughs ("naive" Sachs- Wolfe ef- 
fect) back then. Free-streaming along our (linearly perturbed) past light cone leaves 
the pattern largely unaffected, except that temporal evolution in the gravitational 
potential wells as the photons propagate through them leaves a further AT im- 
print, called the integrated Sachs- Wolfe effect. Intense theoretical work over three 
decades has put accurate calculations of this linear cosmological radiative transfer 
on a firm footing, and there is a speedy, publicly available and widely used code for 
evaluation of anisotropies in a variety of cosmological scenarios, "CMBfast" [9], in- 
cluding the latest hydrogen/helium recombination evaluations, and with extensions 
to more cosmological models added by a variety of researchers. 

Of course there are a number of nonlinear effects that are also present in the maps. 
These secondary anisotropies include weak-lensing by intervening mass, Thompson- 
scattering by the nonlinear flowing gas once it became "reionized" at z ~ 20, the 
thermal and kinematic SZ effects, and the red-shifted emission from dusty galaxies. 
They all leave non-Gaussian imprints on the CMB sky. 

The Future, beyond 2000: We are only at the beginning of the high precision 
CMB era. HEMT-based interferometers are already in place taking data: the 
VSA (Very Small Array) in Tenerife, the CBI (Cosmic Background Imager) in 
Chile, DASI (Degree Angular Scale Interferometer) at the South Pole, where the 
bolometer-based single dish ACBAR experiment will operate this year. Other 
LDBs will be flying within the next few years: Arkeops, Tophat, Beast/Boost; 
and in 2001, Boomerang will fly again, this time concentrating on polarization. 
As well, MAXIMA will fly as the polarization-targeting MAXIPOL. In April 2001, 
NASA will launch the all-sky HEMT-based MAP satellite, with 12' resolution. 
Further downstream, in 2007, ESA will launch the bolometer+HEMT-based Planck 
satellite, with 5' resolution. 

Secondary anisotropies are also being targeted with new instruments. SZ 
anisotropies have been probed by single dishes, the OVRO and BIMA mm arrays, 
and the Ryle interferometer. A number of planned HEMT-based interferometers 



being built are more ambitious: AMI (Britain), the JCA (Chicago), AMIBA (Tai- 
wan), MINT (Princeton). As well other kinds of bolometer-based experiments will 
be used to probe the SZ effect, including the CSO (Caltech submm observatory) 
with BOLOCAM on Mauna Kea, ACBAR at the South Pole, the LMT (large mm 
telescope) in Mexico, and the LDB BLAST. Anisotropies from dust emission from 
high redshift galaxies are being targeted by the JCMT with the SCUBA bolome- 
ter array, the OVRO mm interferometer, the CSO, the SMA (submm array) on 
Mauna Kea, the LMT, the ambitious US/ESO ALMA mm array in Chile, the LDB 
BLAST, and ESA's FIRST satellite. About 50% of the submm background has so 
far been identified with sources that SCUBA has found. 

The CMB Analysis Pipeline: Analyzing Boomerang and other experiments 
involves a pipeline that takes (1) the timestream in each of the bolometer channels 
coming from the balloon plus information on where it is pointing and turns it into 
(2) spatial maps for each frequency characterized by average temperature fluctua- 
tion values in each pixel (Fig. 1) and a pixel-pixel correlation matrix characterizing 
the noise, from which various statistical quantities are derived, in particular (3) 
the temperature power spectrum as a function of multipole (Fig. 1), grouped into 
bands, and two band-band error matrices which together determine the full like- 
lihood distribution of the bandpowers [10]. Fundamental to the first step is the 
extraction of the sky signal from the noise, using the only information we have, the 
pointing matrix mapping a bit in time onto a pixel position on the sky. 

There is generally another step in between (2) and (3), namely separating the 
multifrequency spatial maps into the physical components on the sky: the primary 
CMB, the thermal and kinematic Sunyaev-Zeldovich effects, the dust, synchrotron 
and bremsstrahlung Galactic signals, the extragalactic radio and submillimetre 
sources. The strong agreement among the Boomerang maps indicates that to first 
order we can ignore this step, but it has to be taken into account as the precision 
increases. The Fig. 1 map is consistent with a Gaussian distribution, thus fully 
characterized by just the power spectrum. Higher order (concentration) statistics 
(3,4-point functions, etc.) tell us of non-Gaussian aspects, necessarily expected 
from the Galactic foreground and extragalactic source signals, but possible even in 
the early Universe fluctuations. For example, though non-Gaussianity occurs only 
in the more baroque inflation models of quantum noise, it is a necessary outcome 
of defect-driven models of structure formation. (Peaks compatible with Fig. 1 do 
not appear in non-baroque defect models, which now appear unlikely.) Though 
great strides have been made in the analysis of Boomerang and Maxima, there is 
intense needed effort worldwide now to develop new fast algorithms to deal with 
the looming megapixel datasets of LDBs and the satellites [1 1] . 

COSMIC PARAMETER ESTIMATION 

Parameters of Structure Formation: For this paper, we adopt a re- 
stricted set of 8 cosmological parameters, augmenting the basic 7 used in [7,8], 



£lk, uJb, oJ c dm, n s, t~c cr 8 }, by one. The vacuum or dark energy encoded in 
the cosmo logical constant £7a is reinterpreted as the energy in a scalar 
field Q which dominates at late times, which, unlike A, could have complex dy- 
namics associated with it. Q is now often termed a quintessence field - see 
http://feynman.princeton.edu/ steinh/ "Quintessence? - an overview" for a ped- 
agogical introduction. One popular phenomenology is to add one more parame- 
ter, wq = Pq/pq, where pq and pq are the pressure and density of the Q-field, 
related to its kinetic and potential energy by p Q = Q 2 /2 + (VQ) 2 /2 + V{Q), 
Pq — Q 2 /2 ~ (VQ) 2 /6 — V(Q). Thus wq = —1 for the cosmological constant. 
Spatial fluctuations of Q are expected to leave a direct imprint on the CMB for 
small £, typically smaller than Boomerang or Maxima are sensitive to. We ignore 
this complication here. As well, as long as wq is not exactly — 1, it will vary with 
time, but the data will have to improve for there to be sensitivity to this, and for 
now we can just interpret wq as an appropriate time- average of the equation of 
state. The curvature energy Q k = 1 — Q tot also can dominate at late times, as well 
as affecting the geometry. 

We use only 2 parameters to characterize the early universe primordial power 
spectrum of gravitational potential fluctuations $, one giving the overall power 
spectrum amplitude V®(k n ), and one defining the shape, a spectral tilt n s (k n ) = 
1 + din V<$>/d In k, at some (comoving) normalization wavenumber k n . We really 
need another 2, Vcwikn) an d n t (k n ), associated with the gravitational wave com- 
ponent. In inflation, the amplitude ratio is related to rit to lowest order, with 
0(n s — rit) corrections at higher order, e.g., [6]. There are also useful limiting 
cases for the n s — n t relation. However, as one allows the baroqueness of the in- 
flation models to increase, one can entertain essentially any power spectrum (fully 
/c-dependent n s (k) and n t (k)) if one is artful enough in designing inflaton potential 
surfaces. As well, one can have more types of modes present, e.g., scalar isocur- 
vature modes (Vi S (k n ) , rii S (k)) in addition to, or in place of, the scalar curvature 
modes (V$(k n ),n s (k)). However, our philosophy is consider minimal models first, 
then see how progressive relaxation of the constraints on the inflation models, at 
the expense of increasing baroqueness, causes the parameter errors to open up. For 
example, with COBE-DMR and Boomerang, we can probe the GW contribution, 
but the data are not powerful enough to determine much. Planck can in principle 
probe the gravity wave contribution reasonably well. 

We use another 2 parameters to characterize the transport of the radiation 
through the era of photon decoupling, which is sensitive to the physical density 
of the various species of particles present then, Uj = Qjh 2 . We really need 4: 
Ub for the baryons, uj c d m for the cold dark matter, wy m for the hot dark mat- 
ter (massive but light neutrinos), and u er for the relativistic particles present at 
that time (photons, very light neutrinos, and possibly weakly interacting prod- 
ucts of late time particle decays). For simplicity, though, we restrict ourselves to 
the conventional 3 species of relativistic neutrinos plus photons, with u er . there- 
fore fixed by the CMB temperature and the relationship between the neutrino 



and photon temperatures determined by the extra photon entropy accompanying 
e + e~ annihilation. Of particular importance for the pattern of the radiation is 
the (comoving) distance sound can have influenced by recombination (at redshift 

z r = a- 1 -I), r s = 6000/^ Mpc ^(c^ + ^er 1 )- 1 / 2 ^ +uj b a/(4uJ 1 /3))- 1 ^ dy/a, 
where u 7 = 2.46 x 10~ 5 is the photon density, u er = 1.68o; 7 for 3 species of massless 
neutrinos and u m = u hdm + u cdm + u h . 

The angular diameter distance relation, 1Z = {d k smh.(xr/d k ),Xr,d k sin(xr/d k )}, 
where Xr = 6000 Mpc f^(uj m + u}qcl~ 6w Q + uo k a)~ 1 / 2 d^fa is the comoving distance 
to recombination, d k = 3000 1 | — x / 2 Mpc is the curvature scale and the 3 cases are 
for negative, zero and positive mean curvature, adds dependence upon u k , ooq and 
wq as well as on u m . The location of the first acoustic peak Lp k is proportional 
to the ratio of 1Z to r s , hence depends upon through the sound speed as well. 
Thus Lpk defines a functional relationship among these parameters, a degeneracy 
[12] that would be exact except for the integrated Sachs-Wolfe effect, associated 
with the change of $ with time if Qq or Q k is nonzero. (If $ vanishes, the energy 
of photons coming into potential wells is the same as that coming out, and there 
is no net impact of the rippled light cone upon the observed AT.) 

Our 7th parameter is an astrophysical one, the Compton "optical depth" tq from 
a reionization redshift z re h to the present. It lowers Ce by exp(— 2tc) at the high 
ts probed by Boomerang. For typical models of hierarchical structure formation, 
we expect t c ^ 0.2. It is partly degenerate with cr 8 and cannot be determined at 
this precision by CMB data now. 

The LSS also depends upon our parameter set: the most important combination 
is the wavenumber of the horizon when the energy density in relativistic particles 
equals the energy density in nonrelativistic particles: kj^ eq ps 5r _1 h _1 Mpc, where 
T ps VL m h.Vt~^ 2 . Instead of V$(k n ) for the amplitude parameter, we often use Cio 
at £ = 10 for CMB only, and of when LSS is added. When LSS is considered in 
this paper, it refers to constraints on T + (n s — l)/2 and ln<7g that are obtained by 
comparison with the data on galaxy clustering and cluster abundances [7]. 

When we allow for freedom in u er , the abundance of primordial helium, tilts 
of tilts (dn{ S;iSjt y(k n )/d\ri k, ...) for 3 types of perturbations, the parameter count 
would be 17, and many more if we open up full theoretical freedom in spectral 
shapes. However, as we shall of now only 3 or 4 combinations can be 

determined with 10% accuracy with the CMB. Thus choosing 8 is adequate for the 
present; 7 of these are discretely sampled [16], with generous boundaries, though for 
drawing cosmological conclusions we adopt a weak prior probability on the Hubble 
parameter and age: we restrict h to lie in the 0.45 to 0.9 range, and the age to be 
above 10 Gyr. 

The First Peak and fi tot , Qq and wq: For given u rn and u>b, we show the 
lines of constant L Pk oc 1Z/r s in the Vt tot -Vtq plane for wq=— 1 in Fig. 2, and in 
the wq-Qq plane for Q tot =l in Fig. 3, using the formulas given above and in [12]. 
Our current best estimate [13] of Lp^., using all current CMB data, is 212 ± 7, 
obtained by forming exp < \nL Pk >, where the average and variance of lnL Pk are 



determined by integrating over the probability-weighted database described above, 
restricted here to the Tc = part. With just the prior-CMB data the value was 
224 ±25, showing how it has localized. The numbers change a bit depending upon 
exactly what database or functional forms one averages over. The constant Lpk 
lines look rather similar to the contours shown in the right panel, showing that the 
lZ/r s degeneracy plays a large role in determining the contours. The contours hug 
the Qfot — 1 line more closely than the allowed Lp^ band does for the maximum 
probability values of uj m and cu b , because of the shift in the allowed L Pk band as 
u m and uj h vary in this plane. 

Marginalized Estimates of our Basic 8 Parameters: Table 1 shows there 
are strong detections with only the CMB data for fltot, and n s in the minimal 
inflation-based 8 parameter set. The ranges quoted are Bayesian 50% values and 
the errors are 1-sigma, obtained after projecting (marginalizing) over all other pa- 
rameters. With Maxima, u c d m begins to localize, but much more so when LSS 
information is added. Indeed, even with just the COBE-DMR+LSS data, uj cdm 
is already localized. That Qq is not well determined is a manifestation of the 
Q tot -Qq near-degeneracy discussed above, which is broken when LSS is added be- 
cause the CMB-normalized <t 8 is quite different for open cf. pure Q-models. Su- 
pernova at high redshift give complementary information to the CMB, but with 
CMB+LSS (and the inflation-based paradigm) we do not need it: the CMB+SN1 
and CMB+LSS numbers are quite compatible. In our space, the Hubble parame- 
ter, h = (^(^jh 2 )) 1 / 2 , and the age of the Universe, to, are derived functions of the 
iljh 2 : representative values are given in the Table caption. CMB+LSS does not 
currently give a useful constraint on wq, though wq < —0.7 with SN1. 

The Influence of Light Massive Neutrinos: In [13], we considered what 
happens as we let Q mu /Q m , the fraction of the matter in massive neutrinos, vary 
from to 0.3, for Boomerang+Maxima+prior-CMB+LSS when the weak-H+age 
+ Vltot = 1 prior probability is adopted. Until Planck precision, the CMB data by 
itself will not be able to strongly discriminate this ratio. Adding HDM does have 
a strong impact on the CMB-normalized as and the shape of the density power 
spectrum (effective Y parameter), both of which mean that when LSS is included, 
adding some HDM to CDM is strongly preferred in the absence of VLq. However, 
though more (cold+hot) dark matter is preferred at the expense of less dark energy, 
significant Qq is still required [15]. The Uf, and n s likelihood curves are essentially 
independent of fl mu /fl m . 

The Future, Forecasts for Parameter Eigenmodes: We can also forecast 
dramatically improved precision with further analysis of Boomerang and Maxima, 
future LDBs, MAP and Planck. Because there are correlations among the physical 
variables we wish to determine, including a number of near-degeneracies beyond 
that for Q tot -QQ [12], it is useful to disentangle them, by making combinations 
which diagonalize the error correlation matrix, "parameter eigenmodes" [6,12]. For 
this exercise, we will add u) hdm and n t to our parameter mix, but set wq—— 1, 
making 9. (The ratio Pgw(^)/^(^h) is treated as fixed by n t , a reasonably 
accurate inflation theory result.) The forecast for Boomerang based on the 440 sq. 



TABLE 1. Cosmological parameter values and their 1-sigma 
errors are shown, determined after marginalizing over the 
other 6 cosmological and 4 + experimental parameters, for 
B98+Maxima-I+prior-CMB and the weak prior used in [7,8] 
(0.45 < h < 0.9, age > 10 Gyr). The LSS prior was also designed 
to be weak. The detections are clearly very stable if extra " prior" 
probabilities for LSS and SN1 are included. (Indeed, they are sta- 
ble to inclusion of stronger priors — except if the BBN-derived 
0.019 ± 0.002 is imposed [7].) Similar tables for B98+DMR are 
given in [7] and for B98+MAXIMA-I+DMR in [8]. If O tot is var- 
ied, but wq = — 1, parameters derived from our basic 8 come out 
to be: age=13.2 ± 1.3 Gyr, h = 0.70 ± 0.09, Q m = 0.35 ± .06, 
fib = 0.065 ± .02. Restriction to O tot = 1 and wq = — 1 
yields: age=11.6 ± 0.4 Gyr, h = 0.80 ± .04, Q m = 0.31 ± .03, 
fib = 0.05 ± .005; allowing wq to vary yields quite similar results. 
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deg. patch with a single 150 GHz bolometer used in the published data is 3 out 
of 9 linear combinations should be determined to ±0.1 accuracy. This is indeed 
what we get in the full analysis CMB only for Boomerang+DMR. If 4 of the 6 150 
GHz channels are used and the region is doubled in size, we predict 4/9 could be 
determined to ±0.1 accuracy. The Boomerang team is still working on the data 
to realize this promise. And if the optimistic case for all the proposed LDBs is 
assumed, 6/9 parameter combinations could be determined to ±0.1 accuracy, 2/9 
to ±0.01 accuracy. The situation improves for the satellite experiments: for MAP, 
we forecast 6/9 combos to ±0.1 accuracy, 3/9 to ±0.01 accuracy; for Planck, 7/9 
to ±0.1 accuracy, 5/9 to ±0.01 accuracy. While we can expect systematic errors 
to loom as the real arbiter of accuracy, the clear forecast is for a very rosy decade 
of high precision CMB cosmology that we are now fully into. 
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FIGURE 1. The Boomerang 150A GHz bolometer map (one out of 16) is shown in 
the top figure. Of the entire 1800 square degrees covered, only the interior 440 sq. 
degs. (within the central rectangle) were used in our analysis, i.e., in all less than 5% 
of the data. The 124 square degree Maxima-I map, drawn to scale, is also shown. In 
the lower figure, the Ci (defined in terms of CMB temperature anisotropy multipoles, 
(AT/T)^ m , as indicated) grouped in bandpowers for prior-CMB experiments, including 
TOCO and the North American Boomerang test flight, (squares) are contrasted with the 
Boomerang-LDB (crosses) and Maxima-I (triangles) results in the upper panel. The lower 
panel shows the optimally-combined power spectra of Boomerang+Maxima+prior-CMB 
(circles) and contrasts it with that for Boomerang+Maxima+DMR (squares), showing 
that including the prior experiments does not make a large difference to the results. 
Best-fit models for arbitrary Sltot an d for 0,^=1 are shown in both panels. 
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FIGURE 2. The left panels show lines of constant Lp^ in the Q m -f2Q plane (assuming 
wq=— 1, i.e., a cosmological constant) for two choices of {w m , a^}, the top the most proba- 
ble values, the bottom when the current BBN constraint is imposed (lowering u>b increases 
the sound speed, decreasing Lp^, and varying u m also shifts it). The 0.65 < h < 0.75 
(heavier shading, Hq) and 11 < age < 15 (lighter shading, to) ranges and decelera- 
tions qo = 0,-1/3,-1/2 are also noted. The sweeping back of the Lp^ curves into 
the closed models as Qq is lowered shows that even if J7 tot =l, the phase space results 
in a ID projection onto the Vttot axis that would be skewed to Vttot > 1, a situation 
we see in Table 1. The right panels show 1,2,3-sigma likelihood contour shadings for 
Boomerang+Maxima+DMR and the weak-H+age prior probability (top left) and when 
LSS is added (bottom left). The supernova contour shadings are also plotted, and the 
solid contour lines are what you get when you combine the two likelihoods. In the bot- 
tom left panel, "prior-CMB" experiments (including TOCO and Boomerang-NA) , have 
been added as well, although it makes little difference to the result whether they are 
included or not. Note that the contours are near the diagonal Qtot = 1 line, but also 
follow a weighted average of Lp^ ~ 210 lines. This approximate degeneracy implies £Iq 
is poorly constrained for CMB-only, but it is broken when LSS is added, giving a solid 
SNl-independent On "detection". 
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FIGURE 3. In the left panels, lines of constant Lpk in the wq-£Iq quintessence 
plane (with Q tot =l) are shown for the most probable values of {w m ,U(,}, and when 
cob is constrained to the BBN value of 0.02. Lines of constant deceleration parameter 
qo = (f2 m + (1 + 3wq)Qq) /2 and 0.65 < h < 0.75 (heavier shading, Hq) and 11 < age 
< 15 (light shading, to) ranges are also shown. The right panels show 1,2,3-sigma like- 
lihood contours for Boomerang+Maxima+prior-CMB data with the weak-H+age prior 
probability and Qf t=l. Top shows CMB only, bottom CMB+LSS. SN1 contours, pro- 
vided by Saul Perlmutter [14], and the CMB+SN1 combined ones are also shown in the 
top, and the CMB+LSS+SN1 ones in the bottom. Thus, little follows from CMB only, 
£Iq gets localized when LSS is added but not wq, and wq localization with SN1 is mainly 
because of SN1. The 1,2,3-sigma lines for SN1 shown in the upper right are rather similar 
to the constant deceleration parameter lines in the left panels. 



